Theory of trapped polaritons in patterned microcavities 
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We consider the system of a quantum well embedded in a planar semiconductor microcavity with 
a shallow circular mesa patterned on top of the cavity spacer. For this system we develop the linear 
coupling theory of polaritons. We then compute polariton eigenstates and the corresponding optical 
spectrum. The theory predicts the existence of laterally confined polariton states with a discrete 
energy spectrum, as well as continuum states above the finite mesa potential barrier. 
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Polaritons in planar semiconductor microcavities 
(MCs) are quasi-particles ideally characterized by free 
motion in two-dimensions. The most remarkable ex- 
perimental signature of free two-dimensional motion is 
the continuous energy-momentum dispersion curve 
More direct evidence was recently provided by measure- 
ments of the polariton motion along the plane 1,(1, 

that 

showed ballistic propagation over hundreds of /im when 
polariton wave-packets were resonantly excited at large 
momentum. 

The free particle picture, however, does not hold ex- 
actly. The structural disorder in the system breaks the 
translational invariance along the plane and produces 
localized polariton states. The polariton localization 
length in a GaAs-based high-quality sample was esti- 
mated in the range of a few tens of /^m close to the band 
bottom, by measuring the momentum broadening of the 
dispersion curve [3, whereas recent measurements of 
spatially-resolved photoluminescence from a ZnSe-based 
sample ^ directly showed the spatial pattern of local- 
ized polariton states. Presumably, polariton localization 
occurs mainly due to disorder in the photonic structure. 
The very small polariton effective mass, as compared to 
the exciton mass, is indeed expected to produce motional 
narrowing that strongly suppresses the effect of exci- 
ton disorder on the polariton states. In the MC, on the 
other hand, small thickness fluctuations can produce size- 
able local variations of the photon resonant energy, that 
are expected to affect the polariton motion. 

The idea of producing a polariton trap by a local varia- 
tion of the MC thickness was recently put forward 8] for 
engineering a "polariton quantum box" , namely a pho- 
tonic structure able to confine polaritons in the three 
spatial directions. Given the very light polariton mass, 
of the order of 10~^ exciton masses, the trap size re- 
quired for producing a sizeable energy quantization is 
300 times larger than in the case of bare excitons. Po- 
lariton quantum boxes as large as 10 um can thus be 
designed. Compared to quantum dots 9], these pho- 
tonic structures are easier to fabricate, allowing perfect 
control over size, shape and positioning. Addressing a 
single quantum box with optical spectroscopy would also 



become a trivial experimental task. Another advantage 
with respect to extended polaritons is the reduced spec- 
tral linewidth - as is the case for quantum dots [lol 
- due to the fact that a spatially localized state samples 
a homogeneous portion of the disordered MC plane. Fi- 
nally, trapped polaritons still retain their quasi-bosonic 
character, contrarily to electron-hole states in quantum 
dots. Indeed, in order for the Pauli exclusion to break the 
bosonic behaviour, it is required that the exciton density 
approach its saturation value [l^ . Populating a single 
conflned state with a few polaritons still keeps the exci- 
ton density well below this limit. This quite unique situa- 
tion, thanks to the discrete energy spectrum, would allow 
Bose-Einstein condensation of polaritons, in contrast to 
the ideally homogeneous two-dimensional case, for which 
the Hohenberg-Mermin- Wagner theorem strictly for- 
bids the occurrence of this quantum phase transition. 

Let us assume a Ac-thick MC with distributed Bragg 
mirrors (DBRs). Following the scheme proposed by El 
Daif et al. in order to create a polariton trap, the 
thickness of the MC spacer is made slightly larger by an 
amount AL within a limited region of the MC plane that 
we call mesa, as sketched in Fig. This local vari- 

ation is obtained by etching the mesa pattern on top of 
the cavity spacer and then growing the top DBR. To an 
increased cavity thickness corresponds a local decrease 
of the resonant photon- mode energy. As the photon cou- 
ples to the exciton, forming a polariton mode, the mesa 
acts as a polariton trap. We point out that only the 
MC spacer thickness is varied, while the optical thick- 
ness of the DBR layers is assumed everywhere equal to 
Ac/4. Therefore, within the mesa the DBRs are slightly 
detuned with respect to the spacer thickness. This makes 
fabrication easier and implies a less pronounced depen- 
dence of the local mode energy on AL. In fact, the 
lowest MC resonant fequency is obtained as the solu- 
tion of (nc/c)[(a; - u;c)Lc + {lo - uJm)LDBR] = 0, where 
Lc is the local cavity thickness, lOc and ujm are respec- 
tively the resonant frequency of the cavity spacer and 
of the DBRs, and Lubr — (Ac/2)nirt2/['T.i(»T.i — ^2)] is 
the field penetration length in the DBRs, expressed in 
terms of the two DBR refraction indices ni > n2- By 
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fects at the mesa contour, Maxwell equations give 
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FIG. 1: (a) Sketched cross section of the mesa structure. The 
various lengths and the number of layers are not in scale, 
(b) Computed energy eigenvalues for the photon modes of a 
circular mesa {D — 8.6 /im, hAcuc = —9 meV). (c) Intensity 
plot of the photon spectral density (log scale, 4 decades from 
blue to red). 



replacing Lc = Xc + AL, — (c/nc)27r/Ac, and uJc — 
(c/nc)27r/(Ac + AL), we obtain Aujc ~ —ujmAL/{Xc + 
Ldbr)- For ideal DBRs tuned to the spacer frequency, 
instead, Auc = —LUmAL/Xc- This weaker dependence on 
AL has thus the further advantage of allowing finer con- 
trol over the energy offset of shallow mesas, for which a 
AL of only a few nm must be obtained in the fabrication 
process. For typical GaAs/AlAs DBRs, Ldbr — 3Ac. 
Given AL = 6 nm in Ref. 0, the previous expression 
predicts a mesa energy offset of hAujc = —9 meV, in 
agreement with the optical characterization of the sam- 
ple. 

For a shallow mesa of lateral extension larger than the 
optical wavelength, we can safely assume that the electro- 
magnetic modes at in-plane position p are locally equiv- 
alent to those of a planar microcavity: 

E(r) -E(p)exp(z/c,(p)z), (1) 

where fcz(p) is piecewise constant. Neglecting border ef- 



V^E(p)-f — eo-fc,^(p) E(p)=0 



(2) 



where Cq is the spacer dielectric constant. Outside the 
mesa, the MC resonance is = 2tt/Xc- Inside the 
mesa we can relate the offset Ak^ to the energy off- 
set Aujc, as Akz = y^Aujc/c. Eq. (O is then solved 
in cylindrical coordinates, assuming a circular mesa of 
diameter D. The eigenmodes are therefore expressed 
as E(p) = Unm{p) exp{im(j>) , where n = 0,1,... and 
m = —n,...,n are the radial and angular mode num- 
bers respectively. The mode eigenenergies and the corre- 
sponding energy-momentum spectral function are plot- 
ted in Fig. n (b) and (c) . A discrete energy spectrum 
appears, whose modes show a flat extended signature 
in momentum space, reflecting their spatial confinement. 
At energy above the 9 meV barrier, we find a continuum 
of states whose energy-momentum signature practically 
coincides to that of a planar MC. This result already sug- 
gests that the structure is able to confine photons in the 
three spatial directions. 

We similarly express the exciton center-of-mass wave 
function in terms of Bessel functions of the first kind, 
V'nm(p) = NnmJmiKnmp) ex.p{im4>), whcrc Knm are the 
exciton eigen-momenta and Nnm is a normalization con- 
stant. These are the modes of a free particle, as the 
exciton motion is not affected by the mesa structure. 
By introducing Bose operators Anm and Bnm for pho- 
ton and exciton modes respectively, the linear exciton- 
photon Hamiltonian can be finally expressed in second 
quantization as 



exc 
nm 



^fit B 



nm-^n'm + b.C. 



(3) 



where ujm'n and uimn^'^ are the eigenenergies of the photon 
and of the (free) exciton modes. As required by symme- 
try, the angular number m is conserved in the coupling. 
The energies are expressed in terms of the Rabi 

splitting of the planar cavity hQn and of exciton-photon 
overlap integrals 



P)- (4) 



For a planar geometry, momentum conservation im- 
plies a one-to-one coupling between exciton and photon 
modes. Here, on the contrary, no selection rule on the 
radial quantum number n exists. For the numerical so- 
lution we therefore choose to retain only a finite number 
of cavity modes Nc and exciton modes for each value 
of m. The resulting polariton modes are obtained by nu- 
merical diagonalization of the {Nc -t- N^) x [Nc + N^) 
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FIG. 2: (a) Computed polariton eigenvalues for D = 3.5 /xm. 
(b) Intensity plot of the polariton energy-momentum spectral 
function for D = 3.5 (log scale, 4 decades from blue to 
red), (c) and (d): same as (a) and (b) for D — 8.6 ^m. 
The dashed lines indicate the bare polariton dispersion of the 
planar MC corresponding to the mesa barrier. 



matrix obtained from the Hamiltonian (PJ. Polariton 
eigenvalues obtained for Z? = 3.5 /im and D ~ 8.6 /im 
are plotted in Fig. |21 (a) and (c) respectively. A detun- 
ing of -1-6.5 meV for the planar cavity mode with respect 
to the exciton was assumed, in order to bring the low- 
est confined photon modes in resonance with the bare 
exciton. The result in Fig. |21 (a) and (c) brings clear 
evidence of a discrete energy spectrum, followed by a 
continuous spectrum at higher energy, both for the lower 
and for the upper polariton. For the smaller mesa, only 
three lower polariton and four upper polariton states are 
present in the discrete spectrum, with energy spacings in 
the order of a few meV. The spectra recently measured 
by El Daif et al. on mesas of comparable size show spec- 
tral features in agreement with the present model. The 
polariton operators obtained from the diagonalization of 
© are expressed as P„m = {X^^B 
Each polariton mode exhibits an angular emission pat- 
tern according to its photon component in momentum 



space, defined as /„„i(k) = |(k|P^„|0)p, that is easily 
computed from the model. By assuming for each mode 
a lorentzian energy spectrum, we can finally compute 
energy-momentum spectral function, as shown in Fig. 
EJb) and (d). In the discrete part of the spectrum, polari- 
ton modes present a flat, broad energy-momentum signa- 
ture, which corresponds to the Fourier transform of spa- 
tially confined states (not shown). The continuous part of 
the spectrum, on the other hand, simply corresponds to 
the dispersion of free two-dimensional polaritons. These 
are the scattering states above the finite energy barriers 
of the potential formed by the mesa. Correspondingly, 
the energy-momentum dispersion is well defined, with 
a negligible broadening in k-space. We point out that 
the discrete modes for D = 8.6 /im are less extended in 
momentum space than those for D = 3.5 /im, forming 
a quasi-continuum that mimicks the energy-momentum 
dispersion pattern of 2-D polaritons. For diameters larger 
than 20 /im, the simulations show a spectrum practically 
identical to that of 2-D polaritons in a MC of thickness 
Ac + AL. 



In conclusion, we have developed a theoretical model 
of polaritons in patterned mesa structures. The model 
predicts a discrete energy spectrum of confined polariton 
states both for the upper and the lower branch, each co- 
existing with a continuous spectrum of extended states 
above the trap barrier. The confined states are extended 
in momentum space - a feature that should be detected 
in angle-resolved spectroscopy. Polariton quantum boxes 
are interesting in several respects. The quantum confine- 
ment is produced in a way that is not as invasive as in 
micropillars [l4| . which seem to considerably reduce the 
quality factor of the photonic structure, presumably by 
opening lateral emission channels. Confined polaritons 
represent a quite unique system of quantum-confined in- 
teracting bosons. They can undergo parametric scatter- 
ing processes - due to the mutual Coulomb interaction 
- that have recently been exploited for producing non- 
classical states of polaritons In quantum boxes, 
momentum conservation no longer holds and parametric 
processes between quantum confined states become al- 
lowed. This might be exploited for designing a quantum 
optical device. Recently, rigurous proof was given, that 
a discrete energy spectrum allows Bose-Einstein conden- 
sation (1^], otherwise impossible in two-dimensions [l3| . 
Polariton quantum boxes, could therefore be the ideal 
system for observing the long sought Bose-Einstein con- 
densation of polaritons. 
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